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ABSTRACT 

Photoelectric emission from dust plays an important role in grain charging 
and gas heating. To date, detailed models of these processes have focused pri- 
marily on grains exposed to soft radiation fields. We provide new estimates of 
the photoelectric yield for neutral and charged carbonaceous and silicate grains, 
for photon energies exceeding 20 eV . We include the ejection of electrons from 
both the band structure of the material and the inner shells of the constituent 
atoms, as well as Auger and secondary electron emission. We apply the model to 
estimate gas heating rates in planetary nebulae and grain charges in the outflows 
of broad absorption line quasars. For these applications, secondary emission can 
be neglected; the combined effect of inner shell and Auger emission is small, 
though not always negligible. Finally, we investigate the survivability of dust 
entrained in quasar outflows. The lack of nuclear reddening in broad absorption 
line quasars may be explained by sputtering of grains in the outflows. 

Subject headings: ISM: dust 

1. Introduction 

Cosmic grains are subjected to electrical charging processes and the resulting non-zero 
charge can have important astrophysical consequences. In Galactic H I regions, the grain 
charge is primarily determined by a balance between starlight-induced photoelectric emission 
and the collisional capture of electrons from the gas. Furthermore, in these regions the gas 
heating is dominated by photoelectric emission from dust. Thus, there have been numerous 
detailed studies of UV-induced photoelectric emission (Spitzer 1948; Watson 1972; de Jong 
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1977; Draine 1978; Tielens & Hollenbach 1985; Bakes & Tielens 1994; Weingartner & Draine 
2001b). 

Photoelectric emission from grains exposed to extreme ultraviolet (EUV) and X-ray 
radiation has received much less attention. Dwek & Smith (1996) modelled, in detail, the 
dust and gas heating for neutral grains exposed to high-energy radiation, but did not address 
the grain charging. 

Charging by high-energy photons can be very important for grains immersed in hot 
plasma and exposed to a hard radiation field. The grain potential is limited by the highest 
photon energy present in the incident radiation; as this increases, the grain can reach higher 
potentials. The efficiency of the gas heating decreases as the grain charge increases, since 
the photoelectrons have to climb out of the potential well. 

Hard radiation fields can be found, e.g., in planetary nebulae and active galactic nuclei 
(AGNs). If the gas contains dust, then a fraction of the incident radiation will be processed 
into infrared (IR) radiation, i.e., thermal dust emission. Dwek & Smith (1996) were moti- 
vated by the prospect of learning about these environments through the analysis of the IR 
emission. Ferland et al. (2002) stressed that the X-ray and IR spectra in AGNs may be 
correlated as a result of the interaction between hot grains and gas. Detailed modelling of 
these regions will only be possible once the grain charging is properly modelled, since the 
grain and gas heating rates depend on the charging. 

In a previous paper (Weingartner & Draine 2001b, hereafter WD01), we modelled in 
detail the photoelectric emission from grains exposed to UV radiation; this is briefly sum- 
marized in §2. In §§3 through 8, we extend the WD01 model to include EUV and X-ray 
radiation. Applications to planetary nebulae and quasar outflows are described in §9 and 
§10, respectively. 

2. Photoelectric Emission Induced by Low-Energy Photons 

As a starting point for this study, we will adopt the WD01 photoelectric emission model. 
Here, we briefly discuss some of the important features of this model and mention two minor 
modifications. 

The photolectric emission rate depends on the absorption cross section and the photo- 
electric yield Y (i.e., the probability that an electron will be emitted following the absorption 
of a photon), both of which depend on the photon energy hu. 

WD01 assumed spherical grains, for which the absorption cross section is Qabs^a 2 , where 
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a is the grain radius and Qabs is the absorption efficiency factor. They used Mie theory to 
compute Qabs and adopted dielectric functions from Li & Draine (2001) and Weingartner & 
Draine (2001). Here, we will use Mie theory to compute Q ahs when x = 2na/\ < 2 x 10 4 
(A is the wavelength of the radiation) and anomalous diffraction theory when x > 2 x 10 4 
(Draine 2003). We also will use somewhat newer dielectric functions (Draine 2003). 

The WD01 prescription for determining the threshold photon energy for photoelectric 
emission, hu pet , is discussed in their §§2.2 and 2.3.1 and equations 2 through 7. Here, we will 
adopt a different expression for the minimum energy that an electron must have to escape 
a negatively-charged grain, E m \ n . Instead of equation 7 from WD01, we adopt equation 1 
from van Hoof et al. (2004). 1 Thus, his pet for photoelectric emission from the band structure 
is taken to be 

lPv(Z,a) ,Z>-1 



hu pet (Z,a) = < (1) 
[IP v (Z,a) + E min {Z,a) , Z < -1 

with the "valence band ionization potential" 

IP V (Z, a) = W+(z + l) e - + (Z + 2) (2) 
\ 2 J a a a 

and 

JO , Z > -1 

Em[n= \e v {v = \z + \\) [i-o.3(a/ioA)-°- 45 |z + i|-°- 26 ] ,z<-i ; (3) 

W is the work function and 9 V is defined in equation 2.4 of Draine & Sutin (1987). We retain 
the WD01 estimates of W = 4.4 eV (8 eV ) for carbonaceous (silicate) grains. 

WD01 adopted the following expression for the photoelectric yield of a grain with charge 
Ze (their equation 12; e is the proton charge): 

Y(hv, Z, a) = y 2 (hu, Z, a) x min [y (Q)y 1 (a, hu), 1] , (4) 

where 

Q= (hu- Hot + (Z + l)e 2 /a ,Z>0 

{his- Ket , Z < U 

(WD01 equation 9). For a bulk solid, Y{hv) = yo(® = hu — W). WD01 present approxima- 
tions, derived from laboratory measurements, for yo(Q) for carbonaceous and silicate grains 
(their equations 16 and 17, respectively). 



1 Since van Hoof et al. (2004) define -E m i n as the negative of our E mln , the minus sign in their equation 1 
is absent here. 
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We adopt the size-dependent yield enhancement factor yi used by WD01: 

, u \ ( « 2 -2a + 2-2exp(-a) 

H = U J ^-2/? + 2-2exp(-/?) (6) 

where (5 = a/l a and a = a/l a + a/l e . It depends on the photon attenuation length, l a (the 
e-folding length for the decrease of radiation intensity as it propagates into the material) and 
the electron escape length, l e (roughly the distance that the electron travels in the material 
before losing its energy). The electron escape length depends on the energy of the excited 
photoelectron; WD01 adopted l e — 10 A in all cases, which is a reasonable approximation 
for the low-energy electrons excited by visible and UV radiation. 

Finally, the factor y 2 (hu,Z,a) accounts for the fact that not all of the electrons that 
breach the surface barrier have sufficient energy to escape to infinity. WD01 assumed a 
parabolic electron energy distribution: 

fO(jp\ 6(E - Ei ow )(E high - E) 

JE\ h ) = 7T, p ^3 » ^low < E < ^ high , {I) 

l-C'high — -C'lowJ 

where f%{E)dE gives the fraction of attempting electrons with energy (with respect to 
infinity) between E and E + dE. When Z < 0, E\ ow = E min and -Ehigh = -Emm + /i^ — /if pe t; 
when Z > 0, i?i ow = — (Z + l)e 2 /a and E high = hu — hu pct . The fraction of attempting 
electrons that escape to infinity is given by 



i,Jhu 7, n\- S icf hlSh dE fE( E ) = ^hi g h(^high - 3-Eiow)/ (^high - ^low) 3 , Z > 

' ' } ll ,Z<0 



(8) 



3. Photoelectric Emission Induced by High-Energy Photons 

Photoelectric emission induced by high-energy photons differs from that induced by 
low-energy photons in the following ways: 

1. Whereas low-energy photons can only excite electrons from the band structure of the solid, 
high-energy photons can also excite inner shell electrons. Except very near the absorption 
edge (see, e.g., Draine 2003) the absorption cross section from inner shell electrons will be 
essentially identical to that for isolated atoms. Dwek & Smith (1996) find that the transition 
between band-like and atomic-like absorption typically occurs for photon energies around 
50 eV . However, laboratory-measured photoelectron yields are only available for hu up to 
w 20 eV . 

2. A sufficiently energetic photon can excite a photoelectron from atomic shells other than 
the highest occupied shell (e.g., from the Is shell in C if hu > 291 eV). The hole left by 
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the photoelectron can then be filled in a radiationless transition in which a second electron 
fills the vacancy produced by the photoionization, with the excess binding energy going 
into kinetic energy of a third electron, which then leaves the atom (the Auger effect). The 
Auger electron can then possibly leave the grain. This process can give rise to a cascade of 
secondary Auger transitions that fill holes produced in previous Auger transitions. Thus, 
more than one electron can be emitted from a grain following the absorption of a high-energy 
photon. 

3. Following the absorption of a high-energy photon, the photolectrons and Auger electrons 
might have enough energy to excite secondary electrons from the grain. 

In the following sections, we will discuss the emission of primary photoelectrons, Auger 
electrons, and secondary electrons. 



4. Emission of Primary Photoelectrons 
4.1. Bulk Yields 

4-1.1. Physical Model 

In the case of high-energy radiation, we assume atomic-like absorption and estimate the 
bulk yield of photoelectrons ejected from each shell of each atomic constituent of the grain: 
yo (i, s; 6j jS ) is the probability that a primary photoelectron excited from shell s of element % 
escapes the bulk solid following the absorption of a photon of energy hv = Oj iS + Jj iS , where 
Jj S is the ionization energy of shell s of element i. 

To estimate the bulk photoelectron yields, we adopt a semi-infinite slab geometry and 
assume that a photoelectron escapes the solid with probability 0.5 exp(— z/l e ), where z is 
the perpendicular distance between the surface and the point of excitation. The factor 0.5 
accounts for the fact that only half of the photoelectrons will travel towards the surface of 
the solid, if they are emitted isotropically. For simplicity, we ignore refraction and the fact 
that the reflectivity varies with the angle of incidence, 9. For an isotropic incident radiation 
field, 



1 cos 9 



-x — 



la le 



'H / iO~i,sle 



r 1 f°° dx 

y (i,s;<d i)S ) = ni<JiJ a / d9 sin 9 cos 9 / — exp 

JO JO la 

(9) 

x is the distance along an incident ray from the surface and ni<j itS l a is the probability that 
absorption is by shell s of element i. In evaluating o- it J a , the photon energy hv = Q i s + I i s ] 
in evaluating l e , the initial energy of the excited photoelectron E e = hv — I i s = Q i s . To find 



le n / -i , la 



7- In 1 + T- 
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the bulk yield of photoelectrons ejected from the band structure, yo(band, 0band), we sum 
Vo(i,s) over all of the shells that comprise the band, except that we take E e = ©band and 
hv = ©band + W. We take photoionization cross sections from Verner & Yakovlev (1995) 
and Verner et al. (1996), making use of the FORTRAN routine phfit2. 2 

The photon absorption length l a = A/(47iTmm), where A is the wavelength in vacuo 
and m(A) is the complex refractive index, can be approximated by 

c 1 (h = J]^m(h • ( io ) 

i,s 

For consistency with our atomistic approach to photoelectric emission, we will use the ap- 
proximation when hv > 20 eV . 



4-1.2. Electron Escape Length 

WD01 took l e = 10 A for hv < 20 eV, in approximate agreement with experiments on 
C (Martin et al. 1987) and Si02 (McFeely et al. 1990) films. For high initial energies E e , the 
electron escape length is roughly given by (Draine & Salpeter 1979) 

/ \ -0.85 / rp \ 1.5 

Z e (£ e )« 300A( — [^) , 300eV < £ e < 1 MeV , (11) 

where p is the density of the material. Extending equation (11) below 300 eV, we would 
find l e — 10 A at E e = 164 eV for carbonaceous grains (assuming the ideal graphite density 
of p = 2.24 gem -3 ) and at E e = 211 eV for silicates (assuming p = 3.5 gem -3 ). The upturn 
in l e at ~ 200 eV for C is in rough agreement with the results of Martin et al. (1987), who 
considered E e up to 1 keV . Thus, we adopt 

, , v f 10 A , E e < 164eV 

l e (E e ; carbonaceous) = | ^ ?g x 10 _ s A( ^ /eV) i. 5 > ^ > 164eV ( 12 ) 



and 



, , . v , . 10 A , E e < 211 eV /1Q , 

sdicate) = <{ 32? x l0 _ 3 HEe/ey)1 , )Ee>meV • (13) 



2 The subroutine phfit2 was written by D. A. Verner and is available at 
http://www.pa.uky.edu/~verner/fortran.html. 
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4-1.3. Carbonaceous Grains 

We adopt the ideal graphite density of 2.24 gem -3 ; thus, the C atom number density 
nc = 1.12 x 10 23 cm -3 . There are strong correlations between the 2s and 2p shells, with no 
jump in the photoionization cross section at the 2s threshold (Verner et al. 1996). Thus, 
we treat (2s + 2p) as a single shell, with ionization energy equal to that of the 2p shell 
= 11.26eV. (See Table 1 for the ionization energies of relevant elements.) The bulk yield 
computed using equation (9) is displayed as the short-dashed curve in Figure 1. Since the 
(2s + 2p) shell produces the band structure, this high-energy yield curve should connect 
continuously with the low-energy curve from WD01 (their equation 16), which is displayed 
as the long-dashed curve in Figure 1. To enforce this continuity, we adopt the yield computed 
here when hv > 50 eV , the WD01 yield when hv < 20 eV , and interpolate between these 
when 20 eV < hv < 50 eV. The result is displayed as the solid curve in Figure 1. 

4-1.4- Silicate Grains 

For silicate grains, we adopt a stoichiometry approximating MgFeSi0 4 and a density of 
3.5 gem -3 , intermediate between the values for crystalline forsterite (Mg 2 Si04, 3.21 gem -3 ) 
and fayalite (Fe 2 Si04, 4.39 gem -3 ). The atomic number densities are thus «M g = ^Fc = 
n Si = 1.22 x 10 22 cm -3 and n = 4.88 x 10 22 cm -3 . 

As with C (2s + 2p), O (2s + 2p), Si (3s + 3p), and Fe (3d + 4s) do not display jumps 
in the photoionization cross section at the threshold energy of the deeper shell. Thus, we 
treat each of these pairs as a single shell, with ionization energy equal to that of the lower- 
energy shell. The O 2s and 2p, Si 3s and 3p, Mg 3s, and Fe 3d and 4s shells are considered 
to comprise the silicate band structure. As with graphite, we adopt the WD01 yield when 
hv < 20 eV , the yield computed here when hv > 50 eV , and interpolate between these when 
20 eV < hv < 50 eV . The resulting yield is displayed in Figure 2. 

4.2. Size-Dependent Yield 

For photoelectric emission from the band structure of carbonaceous and silicate grains, 
we compute yi(hv, a) and i/2(hv, Z, a) using the WD01 prescription, as modified in §2 above 
(eqs. 6 and 8). We estimate 2/o(@band) as discussed in §§4.1.3 and 4.1.4; 6 ba nd is computed 
using equation (5). Equation (4) is employed to compute the yield, except that the 1 in the 
min function is replaced by the probability Pband that the photon absorption occurs in the 
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band structure, rather than in an inner shell. We estimate 

band 

Pb a nd(his) « l a (hv) 7W7 i>a (/M/) , (14) 

where the sum is over all shells that are taken to comprise the band structure (see §§4.1.3 
and 4.1.4). 

The yield for photoelectric emission from the inner shells is computed in exactly the 
same way, except that the work function is replaced by the appropriate ionization energy 
when computing the photoelectric threshold energy, hu pct (which is needed to evaluate O ijS ; 
see eq. 5). Also, Pband is replaced with P i s = l a ni(7^ s . 

Very small carbonaceous grains are taken to be polycyclic aromatic hydrocarbons (PAHs). 
For these, we expect Y — > 1 as hv — > oo for the electrons associated with the "band struc- 
ture". This applies by construction for the WD01 yields, but does not apply to the yields 
derived here. Thus, when a < 6 A (corresponding to « 100 C atoms; see eq. 1 in WD01), 
we employ y from WD01. When 6 A < a < 13 A, we take 

13 A — a a — 6 A 

l/0,band = 2/0,band(WD01) — | h £/0,band (this WOrk) — j . (15) 



5. Emission of Auger Electrons 

Suppose a primary photoelectron is ejected from shell s of element i. There are a number 
of Auger transitions that can then occur; we will denote them with index j. Adopting the 
model for electron escape from §4.1.1, the bulk yield for the j-ih Auger transition is 



V0,A(h s ,3', ®A;i,s) = PA,i,s,jni<Ji :S l e 



^lnfl + ^ 



la V If 



(16) 



where Pa,i,sj is the average number of electrons ejected from the atom via Auger transition 
j, photon energy hv = Qa-,i,s + h,s is used in evaluating a i>s and l a , and the energy E^^sj 
of the Auger electron is used when evaluating l e . We adopt the values of PA,i,s,j and E^sj 
given in Tables 4.1 and 4.2 of Dwek & Smith (1996). Note that the Auger yield is the average 
number of electrons emitted by the grain via the given Auger transition, rather than the 
probability of electron emission, and may exceed unity. 

For all grain sizes and charges, we take the threshold photon energy equal to Jj iS , even 
though it should really be somewhat less than this, due to the presence of the band structure. 
We use equation (4) to evaluate the Auger electron yield Y A j, except that the 1 in the min 
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function is replaced with YA,i,s,j ;max = Pi,sPA,i,s,j- We also adopt somewhat different values 
of the input energies than for the case of primary photoelectrons. In evaluating yo,A, we 
take 0A ; i,s = hu — I i>s , regardless of grain size and charge. In computing yi(hu,a), we 
take E e = Ea,%,s,3 (when evaluating l e ). We use equation (8) to compute yi{hv,Z,a) with 
-Ehigh = E AyiyS j — (Z + l)e 2 /a; E\ ow = E min when Z < and E\ ow = —(Z + l)e 2 /a when 
Z > 0. For Auger electrons, it can occur that £hi g h < 0; in this case, we set y<i = 0. 



6. Emission of Secondary Electrons 

The secondary electron yield is defined as the average number of secondary electrons 
emitted per absorbed photon; it may exceed unity. This yield can be expressed as a sum of 
partial yields (denoted with index k), with a term for each type of process that can excite a 
secondary electron (namely, primary emission from the band structure and inner shells and 
Auger emission): 

Y scc (hv, Z,a) = ^2 Y sec,k( hu i Z,a) . (17) 
k 

Since secondary electron energies are generally low, only those excited close to the surface 
have significant probability of escape. Thus, we take Y sec>k — if Y k — 0; otherwise, we 
approximate 

y 2 ,k{his,Z,a) e a" 1 + (10 A)- 1 

where E e>k = hu — W for primary electrons excited from the band structure, hv — l i)S for 
primary electrons excited from inner shell (i, s), and EA,i, s ,j for an Auger electron; e is 
the average energy loss by the exciting electron per secondary electron created in the solid. 
Draine & Salpeter (1979) estimate that e = 117 eV for graphite and 155 eV for lunar dust; 
we will adopt these values for carbonaceous and silicate grains. The grain size a is included 
in the last factor to account for the fact that the exciting electron may escape the grain 
before producing secondary electrons. 

For the secondary electron energy distribution, we adopt 

SUE) = aE-HE-E^) (lg) 
1 + E--'(E - 

(Draine & Salpeter 1979) with E 2 = 8eV 2 and normalization factor 



a = <1- 



1 + E~' 2 (.Bhigh — E\ ow Y V2 } . (20) 
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Thus, 



yTk(hu,Z,a) 



a 



1 + E~ 2 El 



low 



-1/2 



1 + E 2 (-Ehigh — E\ ow y 



'1/2 



Z > 
Z < 



(21) 

when -Ehigh > 0; otherwise, y^ffc = 0- F° r secondary electrons excited by a primary photoelec- 
tron, we take E\ ow and -Ehigh to be the same as for the primary; likewise for those excited by 
Auger electrons. Consequently, the threshold photon energies for the emission of secondary 
electrons are equal to those for the electrons that excite them. 

Figure 3 displays the yields of primary, Auger, and secondary electrons for carbonaceous 
grains with a = 0.1/im. The primary yield Y p reaches unity at hv = 10 4 eV, where both 
l a and l e are very large. However, Qabs^p is a decreasing function at high photon energy. 
Figures 4 and 5 display the total yield (primary plus Auger plus secondary) for carbonaceous 
and silicate grains of various sizes. 



7. Total Photoelectric Emission Rate and Grain Charging 

The primary photoelectric emission rate (excluding that due to excess "attached" elec- 
trons on negatively charged grains) is 

J pc = HO 2 I dv—^Q &hs F p; band + Y] ^p;i, S (22) 

(c.f. eq. 25 in WD01); hu max is the highest photon energy in the incident radiation field and 
Y p denotes the yield of primary electrons. For negatively charged grains, the excess electrons 
undergo photodetachment at a rate 



Jpd 

'^pdt 



/^ max cu u . . 

L iv ^' • (23) 



the photodetachment threshold energy and cross section are taken from equations 18 and 19 
in WD01, respectively. The Auger electron emission rate is 

/^ max CU -\ 

dv-^Qabs Y A;i,s,j , (24) 
- is, min/h i^gj 

where /j S)m i n is the lowest ionization energy among those of the inner shells. Of course, for all 
shells the yields are zero when hv < hu pet . The emission rate of secondary electrons (excited 
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by primary photoelectrons and Auger electrons) is 



I^max CUy X : 

■Jsec, pe = ltd I dv-T Qabs / ^sec, k i (25) 

./Kpet(band) nu k 

where the sum over k denotes the sum over all exciting primary and Auger electrons. 

The collisional charging rates due to positive ions (assumed to have charge e) and 
electrons are denoted J ion and J e , respectively, and are computed as described in §§2.4 and 3 
in WD01. Since we will be considering high-temperature gas here, we also include secondary 
emission induced by incident gas-phase electrons. We compute the charging rate J S cc,gas 
using the Draine & Salpeter (1979) model (their eqs. 14 through 16 and A7 through All). 

When each charging event involves the transfer of a single charge quantum ±e, statistical 
equilibrium yields 

f z (Z)[J pc (Z) + J A (Z) + J scc , pc (Z) + J pd (Z) + J ion (Z) + Jscc, g as] = fz(Z + l)J e (Z + l) , (26) 

where fz(Z) is the probability for the grain charge to be Ze. WD01 used this equation 
(without Ja, «/ S ec,pc, and J SC c,gas) to solve for the grain charge distribution. When Auger 
and secondary emission are important, more than one electron can be emitted at a time, 
yielding a slightly more complicated balance equation and a very complicated combinatorics 
problem. To simplify matters, we will continue to use equation (26). We also compute a 
single "equilibrium charge", i.e., the charge for which the net current J pc + J p( j + J a + ^scc, pe + 
Jion + </scc,gas — J e equals zero. For large grains, the equilibrium charge is nearly identical to 
the average over the charge distribution. Even for small grains, the difference is not large 
when the grain is highly charged. In the following sections, we will present results obtained 
by resolving the charge distribution, unless otherwise indicated. 

A rigorous treatment of multi-electron ejection events would increase the value of Z 
when computing potentials and yields for successive escaping electrons. We will ignore this 
complication as well. 



8. Gas Heating 

The photoelectric gas heating rate per grain with charge Ze, due to primary photoelec- 
trons originating in the band structure, is given by 



/•"'max „ u /•£ max (band) 

T' pe (a) = na 2 du-^Q ahs Y p . hand d£/ B (band; E)E , (27) 

(band) > lu J_B min (band) 
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where E min is given by equation (3), E max = hv - hu pet + E min , and f E (E) = f%(E)/y 2 . 
The heating rate due to primary electrons originating in inner shell (i, s), s , is identical, 
except that each quantity in equation (27) is evaluated for (i, s). Similarly for the heating 
rate due to Auger and secondary electrons, except that E max = I i s — [Z + l)e 2 /a for Auger 
emission and secondary emission excited by Auger electrons; ij iS is the ionization energy of 
the shell in which the primary excitation occurs. The heating rate due to photodetachment 
is given in eq. 40 of WD01. 

The gas cooling rate per grain due to recombination of charged particles with the grain, 
Ag r , is computed as described in §5 of WD01. The cooling is reduced when secondary 
electrons are ejected. We make a simple estimate of this reduction by assuming that the 
average energy (with respect to infinity) of escaping secondary electrons is 

/if"' dEf° E (E) 

where (Eq) is the average energy of the incident electrons (eq. 16 in Draine & Salpeter 1979), 
is the grain potential, f% is given by equation (19) with E\ ow = and -Ehigh = (Eq), and 
Ei = when Z < and e<f) when Z > 0. 



9. Planetary Nebulae 

A planetary nebula is exposed to a hard radiation field from the hot central star, and 
photoelectric emission from dust is an important mechanism for heating the gas (Dopita & 
Sutherland 2000; van Hoof et al. 2004). 

In Table 2, we display photoelectric heating rates r pe and collisional cooling rates A 
for a model planetary nebula and several different grain size distributions. We assume that 
the radiation field is a blackbody with a color temperature T c = 2 x 10 5 K and intensity 
G = 2 x 10 3 . (G = M^ d /-uj^ b , where u^ d is the energy density in the radiation field between 
6eV and 13.6 eV and u^[ ah = 5.33 x 10~ 14 ergcm~ 3 is the Habing (1968) estimate of the 
starlight energy density in this range.) This intensity corresponds to a distance of 10 17 cm 
from a star with luminosity 10 4 L & . We adopt a fully ionized H nebula with number density 
n H = 10 3 cm" 3 and temperature T = 10 4 K. 

The dust in a planetary nebula is unlikely to be the same as the average interstellar dust 
modeled by Weingartner & Draine (2001a). We assume that either carbonaceous grains or 
silicate grains are present, but not both. The size distributions for either the carbonaceous 
grains or the silicate grains are taken from the size distribution for the appropriate compo- 



-13- 



nent in various dust models developed by Weingartner & Draine (2001a). These reproduce 
extinction laws with Ry = Ay/E(B — V)= 3.1, 4, or 5.5, for various values of the amount 
be of carbon per interstellar H atom in PAHs, with cases "A" and "B" referring to different 
assumptions made concerning the total amount of material in the interstellar grain model. 
The photoelectric heating contribution T pe / (Gn-n) of the silicate component varies from 3.6 
to 23 x lO^^ergs^" 1 - more than a factor 6 - showing the sensitivity to the actual size 
distribution. The heating and cooling rates in Table 2 should be taken as an indication of 
the general magnitude that r pc and A can have. 

For carbonaceous dust, use of the WD01 yields overestimates the heating rate by a 
factor of 1.1 to 2 and the cooling rate by as much as 28% (for the first entry in Table 2). 
Inner shell and Auger emission are unimportant, and secondary emission contributes only 
at the 5% level or less. For silicate dust, the heating (cooling) rate is overestimated by a 
factor of w 1.9 to 2.5 (1.3 to 1.5) when the WD01 yields are employed. Inner shell and Auger 
emission account for about 10% of the total heating and secondary emission contributes at 
the 5 to 10% level. For both types of dust, the heating (cooling) results differ by at most 
5% (1%) when secondary emission induced by incident gas-phase electrons is ignored. 

Since the adopted photoelectric yields are highly uncertain, we also computed the heat- 
ing and cooling rates with the bulk yields increased and decreased by a factor 2. For car- 
bonaceous grains, the maximal changes in the heating (cooling) rates were ~ 30% (20%). 
For silicate grains, the revised yields resulted in changes of factors ~ 1.5 and 0.6 for the 
heating rate; the cooling rate was somewhat less sensitive to the changes in the yield. 

10. Quasar Outflows 

X-ray observations of broad absorption line quasars (BALQSOs) reveal large columns 
of absorbing material (e.g., Gallagher et al. 2002). Optical observations, on the other hand, 
indicate that the nucleus is not substantially reddened (e.g., Reichard et al. 2003). Thus, 
it appears that the outflows either are dust-free or contain only very large (by interstellar 
standards) grains. One possible explanation is that the dust is destroyed by sputtering in 
the hot gas. In this process, energetic gas-phase ions collide with the grain and eject surface 
atoms into the gas. There are two limiting regimes for sputtering: "thermal" sputtering, 
where the drift velocity of the grain is negligible compared to the thermal motions of the 
ions, and "drift" sputtering, where the relative velocity of impinging ions is due primarily 
to motion of the grain through the gas. In the presence of a hard radiation field, the 
grains may acquire large, positive electric potentials, suppressing destruction by sputtering 
(Mathews 1967; Laor & Draine 1993). On the other hand, for sufficiently large potentials, the 
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electrostatic stress exceeds the tensile strength of the grain material and the grain explodes. 
In this section, we estimate grain potentials and lifetimes for conditions likely to characterize 
BALQSO outflows. 



The normalization constant a = 3.07eV -1 w(> 13.6 eV) = 2.26 x 10~ 10 ergeV _1 n 7 (> 
13.6 eV), where u(> 13.6 eV) [n 7 (> 13.6 eV)] is the total energy density [photon num- 
ber density] beyond 13.6 eV. We consider gas temperatures T = 10 4 , 10 5 , and 10 6 K and 
ionization parameters U = n 7 (> 13.6eV)/nn = 0.1, 1, 10, and 10 2 . 

Figures 6a and 7a display the grain potential for carbonaceous and silicate grains as a 
function of size. For comparison, Figures 6b and 7b display the potentials computed using the 
WD01 photoelectric yields; in this case, we only include photoelectric emission from the band 
structure, photodetachment, and collisional charging (including secondary emission induced 
by colliding gas-phase electrons). The use of our new, more realistic yields for photoelectric 
emission from the band structure substantially reduces the computed potential, despite the 
additional charging due to the emission of inner shell and Auger electrons (which account 
for at most m 30% of the charging) and the emission of photo-induced secondaries (which 
account for at most 2% of the charging). Secondary emission induced by gas-phase electrons 
contributes at most at the 5, 10, and 25% levels for U = 10 2 , 10, and 1, respectively. For 
U = 0.1 and T = 10 6 K, this process dominates the grain charging. The "notches" at the 
far left in Figure 6a arise because we adopt higher yields for photoelectric emission from the 
band structure of carbonaceous grains with a < 13 A (eq. 15). To test how sensitive the 
results are to uncertainties in the photoelectric yield, we computed the potentials for the 
case that U — 10 and T = 10 5 K but with the adopted bulk yields multiplied by a factor of 
2 or 0.5. The largest changes in the potential were by factors 1.3 and 0.75. 



10.1. Grain Charging 



We adopt the unobscured quasar spectrum from Sazonov et al. (2004), given by 




(29) 
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10.2. Grain Destruction 

10.2.1. Coulomb Explosions 

The electrostatic stress on a spherical grain is S — ((f) /a) 2 /An; if S exceeds the tensile 
strength S max: then the grain fragments (Draine & Salpeter 1979; Waxman & Draine 2000; 
Fruchter et al. 2001). Real grains are not perfectly spherical and the structure of grain 
material is poorly known; thus S^ax is highly uncertain. The maximum grain potential for 
which fragmentation will not occur is 

^ = 106xlo3v (io4^) 1/2 fe) • (30) 

If S max ~ 10 11 dyncm -2 , as measured for ideal materials, then ion field emission, rather than 
Coulomb explosions, may limit the positive grain charge (Waxman & Draine 2000). 

Figures 6a and 7a show the locus for Coulomb explosions for two possible values S max 
of the tensile strength. The results in Figures 6a and 7a show that if S^ax ~ 10 10 dyncm -2 , 
then for ionization parameter U > 10, Coulomb explosions result in a minimum grain radius 
of several x 10 -3 /zm. However, electrostatic stresses reach 10 11 dyn cm -2 only for the smallest 
grains and very high values of U . 



10.2.2. Thermal Sputtering 

The cross section for a collision between an ion with charge zie and energy Ei and a 
grain with radius a and potential is 

2 / 1 - Ziefi/Ei , Ei > zte<j> f . 

a = n 'a x < _ , . (oil 

\0 ,Ei<Zie<i> V ; 

Note that we have neglected the polarization of the grain by the ion. Integrating over a 
Maxwell speed distribution and summing over the various ions in the gas yields the sputtering 
rate for the case that > 0: 

J_da = 2^ _ 3/2 £ ^ 1/2 r» ^ ^ eM _ Ei/kT) L_^f\ Y s i (Ei—Zit<l>) , 

mat p n H Jzm V Ei J 

(32) 

where p is the mass density of the grain material, pm p is the average mass of the grain 
consituent atoms, m 8 is the mass of ion i, rij is the number density of ion i, and Y s ^ is the 
sputtering yield for ion i (i.e., the probability that an atom from the grain is ejected following 
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a collision with a gas-phase ion). For the case that all of the ions are singly ionized, 

(33) 



-L* = exp(-e^T)± * 
nn etc riH etc 



</>=0 



where 



For the Sazonov et al. (2004) spectrum (eq. 29), 

„ 1 (>13.6eV) = 2.00xl0'c m -»( I5If | 5i=T )(^)" 2 , (34) 

where L is the luminosity of the quasar between 1 eV and 10 keV and r is the distance to 
the central source. Combining equations (33) and (34) yields the following expression for 
the lifetime of a grain with radius a against thermal sputtering: 

Ta = 5 .oo yr (^L-) (£) ( m L V' f-LV -QY 1 (expOT)a 

y \0.1fimJ \10j V 10 46 erg s" 1 J \pcj V-10- 4 cm 3 Ayr~V 

(35) 

1 r 

(exp(e<j)/kT)) a = - da' exp[e<j)(a')/kT] . (36) 

U mln 

We take a m i n = 5 A and use the equilibrium charge for each grain size, rather than summing 
over the charge distribution. Tielens et al (1994) reexamined thermal sputtering as a function 
of gas temperature, obtaining rates very similar to those found by Draine & Salpeter (1979). 
We adopt the approximate expression from Tielens et al. (1994) for the thermal sputtering 
rate for neutral grains (their eq. 4.21 and Table 4); these were derived assuming that H, He, 
C, N, and O are all singly ionized and have solar abundances. 

Figure 8 displays the ratio (exp(e0//cT)) a of r ts to the value it would take if the grains 
were uncharged. Figure 9 presents thermal sputtering lifetimes for grains with a = 0.3/iin 
when L = 10 46 ergs _1 and r = 3pc. The horizontal long-dashed line is the outflow timescale 
r flow = r /v ~ 100yr(r/3pc) if the flow speed v ~ 3 x 10 4 kms~ 1 . The grains survive only 
for quite large values of the ionization parameter or low values of the gas temperature. 
Furthermore, when U > 10, grains with a ~ 0.3;um experience supersonic drift relative to 
the gas, increasing the destruction rate. Thus, the lack of nuclear reddening in BALQSOs 
does not appear to require that the outflow originates in a dust-free region (e.g., within the 
sublimation radius). However, better observational constraints on the location and physical 
conditions in the outflow are needed before a firm conclusion can be reached. 



-17- 



11. Summary 

In this paper, we have extended the WD01 photoelectric emission model to higher 
photon energies, enabling the treatment of dust exposed to hard radiation fields. We provide 
revised yields for photoelectric emission from the band structure when hv > 20 eV and 
treat the emission of Auger electrons and secondary electrons excited by primary and Auger 
electrons. 

For HII regions, ionized by OB stars, the radiation field is soft enough that the WD01 
model can be employed without modification; the new treatment yields a < 10% difference 
in the gas heating rate. 

For planetary nebulae, the gas heating rates found with the extended model are lower, 
by a factor 2 or less, than those estimated with the WD01 model, since the revised yields for 
photoelectric emission from the band structure are lower. Inner shell, Auger, and secondary 
emission are not important. For quasar outflows, the new treatment yields significantly lower 
grain potentials. In this case, combined inner shell and Auger emission contribute at as high 
as the 30% level. This is less than the uncertainty associated with photoelectric emission 
from only the band structure. Secondary emission is unimportant. It appears that thermal 
sputtering can destroy dust in BALQSO outflows in less than a dynamical time, for a wide 
range of plausible physical conditions in the outflow. However, better constraints on the 
outflow conditions are needed to definitively settle the issue of dust survivability. 

In environments with harder radiation fields (e.g. supernova remnants and the central 
regions of galaxy clusters), inner shell, Auger, and secondary emission may be more impor- 
tant relative to photoelectric emission from the band structure. However, in these cases the 
gas is so hot that secondary emission induced by gas-phase electrons and ions incident on 
the grain is expected to dominate over photoelectric emission in grain charging (Draine & 
Salpeter 1979). 

Support for this work, part of the Spitzer Space Telescope Theoretical Research Pro- 
gram, was provided by NASA through a contract issued by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. BTD was supported in part 
by NSF grant AST-0406883. We are grateful to Norm Murray, Gary Ferland, and Peter van 
Hoof for helpful discussions. 
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10 1 10 2 10 3 10 4 

hi/ (eV) 

Fig. 1. — Bulk yield for carbonaceous dust. Long-dashed curve: WD01 yield; short-dashed: 
yield computed with eq. (9); solid: our adopted yield 
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Fig. 2. — Bulk yield for silicate dust. Long-dashed curve: WD01 yield; short-dashed: yield 
computed with eq. (9); solid: our adopted yield 
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Fig. 3. — The yields of primary, Auger, and secondary electrons for uncharged carbonaceous 
grains with a = 0.1/xm. 
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Fig. 4. — The total yield (primary plus Auger plus secondary) for uncharged carbonaceous 
grains of four different sizes, as labelled. 
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Fig. 5. — Same as Fig. 4, except for uncharged silicate grains. 
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Fig. 6. — (a) Carbonaceous grain potential <fi versus radius a for grains exposed to an 
unobscured quasar spectrum (eq. 29) with ionization parameter U = 0.1, 1, 10, and 10 2 , 
and for gas temperatures T = 10 4 , 10 5 , and 10 6 K, as indicated. For a given value of U, 
the potential is larger for higher T. Coulomb explosions occur to the left of the heavy 
curve labelled "CE(10)" if the tensile strength S max « 10 10 dyncm~ 2 and to the left of the 
curve labelled "CE(ll)" if S max ~ 10 n dyncm~ 2 . (b) Computed grain potentials using the 
photoelectric yields from WD01, for comparison. 




Fig. 7.- 



Same as Fig. 6, except for silicate grains. 
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Fig. 8. — Ratio of the thermal sputtering timescale to its value for uncharged grains, versus 
gas temperature. 
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T (K) 

Fig. 9. — Thermal sputtering lifetime r ts for grains with a = 0.3yum when L = 10 46 erg s _1 and 
r = 3pc, versus gas temperature. The horizontal long-dashed line is the outflow timescale 
r flow , assuming that the outflow speed v ~ 3 x 10 4 kms _1 . 
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Table 1. Ionization Energies' 



Element Is 2s 2p 3s 3p 3d 4s 



C 291 19.39 11.26 

O 538 28.48 13.62 

Mg 1311 94.0 54.9 7.646 

Si 1846 156 106 15.17 8.152 

Fe 7124 857 724 104 66 14.7 7.902 



in eV 
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Table 2. Photoelectric Heating and Recombination Cooling Rates for a Model Planetary 

Nebula 



composition R v 10 5 b c case r p e/(GnH) a A/(Gnn) A 



carbon 


3.1 


0. 


A 


16.2 


1.97 


carbon 


3.1 


6. 


A 


45.9 


5.56 


carbon 


4.0 


0. 


A 


12.8 


1.55 


r*a rhnn 


4 n 


A 


A 


tJl.U 


3 89 


carbon 


5.5 


0. 


A 


9.75 


1.19 


carbon 


5.5 


3. 


A 


23.4 


2.84 


carbon 


4.0 


0. 


B 


16.7 


2.05 


carbon 


4.0 


4. 


B 


32.0 


3.89 


carbon 


5.5 


0. 


B 


13.3 


1.64 


carbon 


5.5 


3. 


B 


22.6 


2.75 


silicate 


3.1 


0. 


A 


15.6 


1.29 


silicate 


3.1 


6. 


A 


22.5 


1.89 


silicate 


4.0 


0. 


A 


9.84 


0.818 


silicate 


4.0 


4. 


A 


9.31 


0.775 


silicate 


5.5 


0. 


A 


3.83 


0.314 


silicate 


5.5 


3. 


A 


3.61 


0.297 


silicate 


4.0 


0. 


B 


10.1 


0.840 


silicate 


4.0 


4. 


B 


9.73 


0.810 


silicate 


5.5 


0. 


B 


3.87 


0.316 


silicate 


5.5 


3. 


B 


3.84 


0.314 



a 10 25 ergs 1 



